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ABSTRACT 

We present Spitzer Space telescope spectroscopic observations of 14 carbon-rich AGB stars 
in the Small Magellanic Cloud. SiC dust is seen in most of the carbon-rich stars but it is weak 
compared to LMC stars. The SiC feature is strong only for stars with significant dust excess, 
opposite to what is observed for Galactic stars. We argue that in the SMC, SiC forms at lower 
temperature than graphite dust, whereas the reverse situation occurs in the Galaxy where SiC 
condenses at higher temperatures and forms first. Dust input into the interstellar medium by 
AGB stars consists mostly of carbonaceous dust, with little SiC or silicate dust. Only the 
two coolest stars show a 30-micron band due to MgS dust. We suggest that this is due to the 
fact that, in the SMC, mass-losing AGB stars generally have low circumstellar (dust) optical 
depth and therefore effective heating of dust by the central star does not allow temperatures 
below the 650 K necessary for MgS to exist as a solid. Gas phase C2H2 bands are stronger 
in the SMC than in the LMC or Galaxy. This is attributed to an increasing C/O ratio at low 
metallicity. We present a colour-colour diagram based on Spitzer IRAC and MIPS colours to 
discriminate between O- and C-rich stars. We show that AGB stars in the SMC become carbon 
stars early in the thermal-pulsing AGB evolution, and remain optically visible for '--^ 6 x 10^ 
yr For the LMC, this lifetime is 3 x 10^ yr The superwind phase traced with Spitzer lasts 
for ~ lO"* yr Spitzer spectra of a K supergiant and a compact HII region are also given. 

Key words: circumstellar matter - infrared: stars — carbon stars — AGB stars — stars: mass 
loss — Magellanic Clouds 



1 INTRODUCTION 

The late stages of the evolution of low and intermediate mass stars 
(hereinafter LIMS) are characterised by intense mass loss. This 
so-called superwind occurs during the Asymptotic Giant Branch 
(AGB) phase. The mass loss leads to the formation of a circumstel- 
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lar envelope composed of gas and dust. The chemical composition 
depends strongly on the C/O abundance ratio. CO is one of the first 
molecules to form and is very stable and unreactive. If C/O > 1 by 
number, O will be trapped in the CO molecules, leading to a chem- 
istry dominated by molecules such as C2, C2H2, HCN and SiC dust 
grains (carbon stars). On the other hand, if C/O < 1 we will observe 
oxygen-rich stars showing SiO, OH, H2O molecules and silicate 
dust. 
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The study of this mass loss is of high astrophysical interest, 
and impacts on the chemical evolution of galaxies. Indeed, the mass 
loss from LIMS contributes to roughly half of all the gas recycled 
by stars (Maeder 1992). Mass loss from LIMS is one of two main 
sources (along with WR stars and supernovae associated with mas- 
sive stars), of carbon in the universe (Dray et al. 2003). LIMS are 
also the main source of heavy s-process elements (e.g. Ba, Pb) and, 
when including post-AGB evolution (e.g. novae), the major stellar 
source of lithium. 

The mass-loss process is, however, not fully understood. The 
mass loss results from a complicated interplay between stellar pro- 
cesses (turbulent convection, pulsation) and circumstellar processes 
(pulsation-driven shocks, radiation pressure) where especially the 
role of the dust composition is disputed (Woitke 2006). The ef- 
fect of metallicity on the mass-loss rates is poorly understood and 
may vary with dust chemistry. Observational studies show no dif- 
ference between peak gas mass-loss rates for LMC and SMC stars 
(e.g. van Loon 2006). Theoretical studies (e.g. Bowen & Willson 
1991) predict that the mass-loss rates in AGB stars depend on the 
metallicity: a lower metallicity leads to a lower dust-to-gas ratio 
and less efficient dust-radiation pressure. The effect may differ be- 
tween mineral (e.g. silicates) and non-mineral dust. In 0-rich stars, 
all dust species are minerals (containing Si, Al for respectively sili- 
cates and corondum), whilst in C-rich stars dust contains both min- 
eral (e.g SiC) and non-mineral components (e.g. soot, amorphous 
carbon). 

The sensitivity reached by the 5pitzer Space Telescope 
(Werner et al. 2004) enables, for the first time, the determination 
of mass-loss rates from stars of different masses all along the AGB 
sequence at the distance of the Magellanic Clouds. The distances 
to these two galaxies are also relatively well-known and the metal 
abundance of stars within can be estimated using age-metallicity re- 
lations. This allows one to measure absolute (dust) mass-loss rates 
for stars with known bolometric magnitudes and metallicities. 

We have therefore conducted a survey of mass-losing stars 
in the Small and Large Magellanic Clouds (hereinafter SMC and 
LMC respectively). The aim of this project is to empirically cali- 
brate the mass-loss rate of AGB stars as a function of mass, lumi- 
nosity and metallicity. Here we present the data from the SMC. The 
data from the LMC are presented in Zijlstra et al. (2006). A study 
of the mass-loss rates will be presented in forthcoming papers. 



2 TARGET SELECTION 

We selected 17 stars in the SMC to obtain a sample of stars of 
lower metallicity than the LMC stars described in Zijlstra et al. 
(2006). The stars were all AGB stars, with preference being given 
to stars in the populous intermediate-age cluster NGC 419, stars 
with past infrared monitoring showing large-amplitude variability 
or stars with some ISO detection (Cioni et al. 2003). 

The observed stars are shown in the Mk vs J— K diagram in 
FiglU Fourteen of the intended targets were observed with Spitzer, 
while peak-up failures occurred for 3 objects. In one case, the peak- 
up occurred on a noise spike so no target was observed, while in 
the other two cases a brighter object in the peak-up image led to 
the detection of an HII region (MB 88, from Murphy & Bessell 
2000) and a red supergiant of spectral type K (PMMR 52, from 
Prevot et al. 1983). The latter object is the single object lying in the 
supergiant region of Figure[T] 

The observed stars in the cluster NGC 419 whose names con- 
tain LE (NGC 419 LE 16, NGC 419 LE 18, NGC 419 LE 27 and 



-15 



-10 



-5 



Supergiants 



Field star 

MSX or IRAS source 
> SMC Spitzer target 
J Cluster star 



FGB tip 




2 4 6 

J-K 



Figure 1. The Mk, J— K diagram for the observed sample of mid-IR 
sources in the SMC (large filled circles). Stars in the cluster NGC 419 
ai'e circled. For orientation, we also show a large sample of mid-infrared 
sources in the LMC from MSX (small open circles) and field stars from a 
small area of the LMC bar (small dots) - see also Zijlstra et al. (2006). AGB 
stars are confined approximately to the region below the line marked "AGB 
limit". Stars above this limit are red supergiants with masses above ~8 Mq , 
or perhaps foreground stars. Distance moduli of 18.54 and 18.93 have been 
assumed for the LMC and SMC, respectively (Keller & Wood 2006). 



NGC 419 LE 35) were discovered by Lloyd Evans (1980) using 
photographic surveys in the V and I bands. The two other stars 
from this cluster, NGC 419 IRl and NGC 419 MIRl were discov- 
ered during a survey with ISOCAM (Tanabe et al. 1997). LEGC 
105 was discovered during an optical survey of the SMC (Lloyd 
Evans et al. 1988). RAW 960 appeared in the Rebeirot et al. cat- 
alogue (1993). IRAS 00554-7351 is described in Whitelock et al. 
(1989). Stars with names beginning with ISO were selected from 
the ISO/MACHO catalogue of variable AGB stars in the SMC 
(Cioni et al. 2003). The variable star GM780 was discovered by 
Glen Moore using UK Schmidt plates but his work is unpublished. 
Literature data is summarised in Table [T] For LEGC 105, Lloyd 
Evans et al. (1988) suspect that the pulsation is irregular but their 
approximate period is confirmed with MACHO. RAW 960 has a 
very short period and is clearly a semi-regular rather than a Mira 
variable. 

All 14 AGB stars that were observed turned out to be carbon 
stars. The dominance of C-stars is at variance with the MSX-based 
classification scheme of Egan et al. (2001) according to which the 
sample stars would be a mixture of carbon-rich and oxygen-rich 
stars. Our LMC sample of AGB stars also showed that the Egan et 
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al. classification scheme is ineffective at separating M- and C-stars 
on the AGB (Zijistra et al. 2006). 

The AGB stars in our SMC sample have JHK colours which 
are reddened compared to stellar colours, but mostly not extremely 
so. The J— K versus Mk diagram is shown in Fig. 1. The SMC sam- 
ple has predominantly lower circumstellar extinction (bluer J— K 
colours) than the LMC sample in Zijistra et al. (2006). The cluster 
stars from NGC 419 are mostly at the blue end of the AGB se- 
quence, but they also include the reddest star in the sample (NGC 
419 MIRl). IRAS 00554-7351 is also very red with J-K= 5. 

Fig.[2]shows the colour-magnitude diagrams for known optical 
carbon stars in the Clouds. Open circles show the Spitzer samples of 
Zijistra et al. (2006) and this paper. The distribution shows that we 
are observing the peak of the luminosity range of the optical carbon 
star distribution. The LMC sample contains a larger fraction of red- 
der objects with higher dust mass-loss rates than the SMC sample. 
The gas mass-loss rates may be more compatible (e.g. Matsuura et 
al. 2006, van Loon 2006). However, this point should be taken into 
account when comparing the samples, as later in this paper. 
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Figure 4. Spectrum of the red object MB 88. This object shows emission 
lines and could be associated with an Hll region. Emission lines from Siv 
(10.51/xm), Nell (12.81fjm), SlII (18.71 and 33.48/im), Sill (34.83/im) and 
Nelll (35.97/jm) are obsei-ved. 



3 OBSERVATIONS 
3.1 Spitzer 

The observations were made with the InfraRed Spectrograph (IRS, 
Houck et al. 2004), on board the Spitzer Space Telescope. We 
used the Short-Low (SL) and Long-Low (LL) modules to cover 
the wavelength range 5-38/im. The SL and LL modules are each 
divided in two spectral segments, together known as SL2, SLl , 
LL2 and LLl; a "bonus" order covering the overlap between the 
two modules is also available. The data reduction is similar to that 
described in Zijistra et al. (2006). The raw spectra were processed 
through the Spitzer pipeline S 12. We replaced the bad pixels by val- 
ues estimated from neighbouring pixels. The sky subtraction was 
done by differencing images aperture by aperture in SL and nod 
by nod in LL. We used the software Spice (Spitzer IRS Custom 
Extractor) to extract the spectra. The flux calibration was made us- 
ing the reference stars HR 6348 (KO III) in SL and HR 6348, HD 
166780 (K4 III) and HD 1735 11 (K5 III) in LL. The spectra were 
individually extracted from the individual images. Both nods in 
both apertures were then joined simultaneously, recalculating the 
errors in the process by comparing the nods. The different nods 
were averaged, using the differences to estimate the errors. The dif- 
ferent spectral segments were combined using scalar multiplication 
to eliminate the discontinuities due to flux lost because of pointing 
errors. The different segments were also trimmed to remove du- 
bious data at their edges. We also retained the bonus order where 
it was valid. We obtained a standard wavelength calibration accu- 
racy of 0.06/xm in SL and 0.15/im in LL. The calibration process 
is detailed in Sloan et al. (2003). 

Even after sky removal, some interstellar emission lines re- 
main on the spectra, mostly observed longward of 30/im. 

Two of the observed objects in NGC 419, LE 35 and LE 27, 
showed large discrepancies between SL and LL. The apertures of 
these two modules are perpendicular to each other on the sky. An 
acquisition error may have led to the star being missed in one aper- 
ture but not the other, or confusion in the cluster may have caused 
spectrum extraction problems. The SL spectra are those of AGB 
stars, with flux densities consistent with the MSX flux. The LL 
spectra are too faint and featureless. We therefore consider the SL 
observations to have been successful and will only discuss their 
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Figure 5. Spectrum of the oxygen-rich star PMMR 52. Typical spectral fea- 
tures of silicates are observed at 9.8^m and 18/im. 



SL spectra. However, mispointings can affect the wavelength cal- 
ibration, flux calibration and even the slope of the spectral energy 
distribution (SED). We can therefore not be as confident of these 
spectra as we are for the remaining stars. 

The spectra of the observed AGB stars, ordered by dust tem- 
peratures (see Sect. 6), are presented in Fig. [3] The molecular bands 
discussed below show that all 14 objects are carbon-rich stars. We 
calculated equivalent magnitudes of our sources for IRAC Spm and 
MIPS 24/im, by convolving these filter band passes with our spec- 
tra. These are listed in Table|2] 

Finally, we present the spectra of the two objects acquired ac- 
cidently following peakup failure. MB 88 (Fig. |4]l is a red object 
with emission lines, a spatial diameter of about 8 arcsec, and with 
several embedded stars. It is almost certainly an HII region. Emis- 
sion lines from SiV (10.5lMm), Nell (12.81/im), SlII (18.71 and 
33.48/im), Sill (34.83^m) and Nelll (35.97^m) are observed. Its 
spectrum is not analysed further here. The spectrum of the red su- 
pergiant PMMR 52 can be found in Fig.|5] It is discussed briefly in 
Sect. [511 
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Figure 2. The Mk, J— K diagram for all known optical carbon stars in the LMC and SMC. Samples are taken from Kontizas et al. (2001) and Rebeirot et 
al. (1993); infrared magnitudes are from 2MASS. Open circles show the Spitzer targets discussed in Zijlstra et al. (2006) and in this paper LMC photometry 
for these is from Groenewegen et al. (2007): where no J is available, we assumed J - K ~ 0.6 + 1.84 X (H - K). The J-Kfa 7 for NGC418MIR1 was 
estimated from its K— L colour 



3.2 Bolometric magnitudes and initial masses 

We determined the bolometric magnitudes of the observed stars 
assuming a SMC distance modulus of 18.93 ± 0.024 (Keller & 
Wood 2006) (~60kpc). JHKL photometry was taken nearly simul- 
tanously with the Spitzer observations. These observations were 
made using the 2.3m Telescope at Siding Spring Observatory (SSO, 
Australia). The filters used were centred at 1.28^m (J), 1.68/im 
(H), 2.22/im (K) and 3.59/im (L). The observations are described 
in Groenewegen et al. (2007). Zero flux was assumed at the fre- 
quencies and 3 x 10^^ Hz to estimate the bolometric magnitudes. 
The flux from the star was then estimated by integrating under the 
resulting SED. The line-of-sight extinction towards the SMC of 
E(B - V) = 0.12 (Keller & Wood 2006) was ignored. The bolo- 
metric magnitudes are presented in Table |2] 

Simultaneous observations in the K-band are important to 
obtain a reliable instantaneous luminosity estimate. However, the 
magnitudes are single-epoch measurements and therefore subject 
to pulsation-induced variability. Earlier (literature) values are listed 
in Table 1. These generally agree to within 0.5 mag of the current 
measurements, consistent with the variability expected for these 
stars. 

The histogram of the distribution is shown in Fig.|6] The SMC 
sample includes the carbon stars observed by Sloan et al. (2006) 
which cover a similar luminosity range. The LMC carbon-star sam- 
ple of Zijlstra et al. (2006) is also shown. There is no obvious dif- 
ference between the LMC and SMC distributions. 



For comparison, in the third panel we show a complete carbon 
star luminosity function for the LMC. This is derived from the car- 
bon star survey of Kontizas et al. (2001), which we cross-correlated 
with 2MASS. The bolometric magnitudes for this sample are de- 
rived from the 2MASS JHK magnitudes, using the bolometric cor- 
rection equation derived by Whitelock et al. (2006): 

BCk = + 0.972 + 2.9292 x {J - K) - 1.1144 x ( J - Kf 
+0.1595 X (J - Kf - 9.5689 10^(J - (1) 

This relation differs by 0.2-0.3 mag from the one used by Costa & 
Frogel (1996) fori < J — K < 2.5. Positional agreement between 
the 2MASS object and the carbon star is required to be better than 
1 arcsec. Stars with non-detections in one or more infrared band are 
not included. 

The bottom panel shows the same for the SMC carbon stars 
taken from Rebeirot et al. (1993). These stars were also cross- 
correlated with 2MASS, and the bolometric magnitudes derived 
from the 2MASS JHK magnitudes. The coordinates in Rebeirot 
et al. are less accurate and we accepted co-identifications out to 2 
arcsec. A larger number of chance superpositions may be expected. 

If we compare this histogram with Fig. 20 in Vassiliadis & 
Wood (1993), for SMC metallicities, this indicates that the stars we 
observed have initial masses in the range 1-4 Mq. 

Six of the stars we observed are located in the cluster NGC 
419. This cluster has an age of 1.2 x 10^ yr, and a metallic- 
ity [Fe/H]=-0.60±0.21 (de Freitas Pacheco et al. 1998). The 
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Figure 3. Spectra of the carbon-rich stars observed in the SMC. The two stars at the bottom have been observed only in SL. All other stars are ordered by dust 
temperature. Blue stars are at the top and red ones at the bottom. Molecular templates are shown at the top (see Sect. 4.3). 
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Figure 6. Top: distribution of bolometric magnitudes of the observed stars 
in the SMC, assuming a distance modulus of 18.93 (Keller & Wood 2006). 
Both stars from this paper and from Sloan et al. (2006) are included. Sec- 
ond: the histogram for the LMC carbon stars observed by Spitzer (Zijlstra 
et al. 2006). Third: Optical carbon stars in the LMC (Kontizas et al. 2001). 
Bottom panel: Optical carbon stars in the SMC (Rebeirot et al. 1993). The 
arrows indicate the location of the tip of the RGB (BeUazzini et al. 2001). 



Table 1. Bolometric magnitudes and periods from the literature. B83: 
Bessell et al. (1983); W91: Westerlund et al. (1991); T97: Tanab^ et al. 
(1997); vL05: van Loon et al. 2005; W89: Whitelock et al (1989); LE88: 
Lloyd Evans et al. (1988); C03: Cioni et al (2003). The distance modulus to 
the SMC is taken as 18.93. 
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Figure 7. The [6.4]— [9.3] versus [16.5]— [21.5] colour-colour diagram of 
the stars in our SMC sample. 



loss is weaker (or the onset of the superwind later) than assumed 
for the evolutionary model calculations (Zijlstra 2004). 



isochrones of Pietrinferni et al. (2006) for this age and for Z = 
0.004, give an initial mass of the thermal-pulsing AGB stars in this 
cluster of 1 . 82 IVIq for standard models and 2.08 for models us- 
ing overshoot. (The difference shows the considerable uncertainty 
in deriving stellar ages.) But their bolometric magnitudes (sec also 
van Loon et al. 2005) are associated with higher-mass progenitors 
on the VassiUadis & Wood (1993) tracks. Thus the stars reach a 
higher luminosity than predicted. This may indicate that the mass 



4 OBSERVED SPECTRAL FEATURES 

The spectra show clear and deep molecular bands. The two main 
dust features are less obvious. We will describe the dust bands first, 
followed by the observed molecular bands. 

4.1 Dust bands 

Most of the carbon stars in our sample show the presence of an 
emission feature around 1 1 .3//m. This feature is attributed to emis- 
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Table 2. Observed SMC targets: names, adopted coordinates, and photometry. J and K are taken from near-simultaneous measurements at SSO. The 8/im 
and 24^m magnitude are calculated from our spectra using the IRAC and MIPS filter profiles; the zeropoints are taken as 63.5 ( 8.0/im) and 7.14 Jy (lA^vn) 
respectively. 
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NGC 419 LE 27 


01082067- 
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sion from SiC. However, in several cases the feature is very weak, 
and in such cases it is hard to judge whether this feature is an 
emission feature or due to molecular absorptions on the blue and 
red side. SiC condenses at high temperatures and is expected to be 
present for all dusty carbon stars. However, its abundance may be 
limited by the abundance of Si. 

NGC419 LE 35 shows no discemable SiC feature. This is, 
however, one of the two sources with acquisition and/or extraction 
problems. GM 780 has a very strong feature, compared to the other 
stars, but the shape is peculiar. The shape is possibly affected by ex- 
traction problems but the strength of the feature is well established. 
This object is discussed below. 

A wide emission feature, centred at ^ 30/.im is observed in 
two stars: NGC 419 MIRl and (weaker) IRAS 00554-7351. This 
feature is common in Galactic AGB and post-AGB stars, and at- 
tributed to emission from MgS (Hony et al. 2002). These two SMC 
stars have the reddest near-infrared colours in our sample, and the 
coolest dust, with dust temperatures of respectively 409 and 589 
K (see section|6]l. In the LMC, this feature is observed only in the 
envelopes of stars with dust temperature lower than 650K (Zijl- 
stra et al. 2006). This can be explained by the fact that MgS grows 
on the surface of pre-existing grains, this process starting around 
600K and being complete around 300K. The two SMC detections 
are consistent with this. 



4.2 Molecular bands 

Most vibrational bands of simple molecules occur at wavelengths 
shortward of 20/im. (HNC has a vibrational fundamental at 20^m.) 
This is the region where we observe clear bands. The main fea- 
tures, observed for most of the stars of our sample, are located at 
~7.5/im and ~13.7^m respectively. Both are associated with ab- 



sorption from C2H2. The 13.7/im band is due to the 1/5 bending 
mode of acetylene and its associated hot and combination bands 
(Cernicharo et al. 1999). This absorption band is observed in almost 
all carbon stars of our sample, but it is not clear in ISO 00549 and 
GM 780. The narrow band is flanked by broader bands, best seen in 
NGC419 MIRl. The double peaked band at 7.5/im observed in all 
of our spectra is associated with the C2H2 vl+i'l P- and R-branch 
transition, as shown by Matsuura et al. (2006). The shape may be 
different in NGC419 LE 27, perhaps with a contribution from CS 
and/or HCN. 

A weak absorption band at 14.3/im is observed in the spec- 
tra of the two reddest stars of our sample (IRAS 00554 and NGC 
419 MIRl). This band has been observed in Galactic carbon stars 
(Aoki et al. 1999) and is attributed to HCN. The fact that it is ob- 
served in the reddest stars of our sample, where the optical depth is 
important, and dust coolest, seems to indicate that this molecule is 
present in the outer layers of the observed envelopes. 

A decrease of the spectral energy distribution is observed near 
the blue edge of all our spectra , around 5 ^m. This is due to several 
molecules: C2 and/or C3 starting from 6/im (blueward), and CO 
absorption very close to the 5 /im. Several targets show a higher 
flux level at the blue edge with a sharp drop: this is not seen in the 
LMC spectra. 

Zijlstra et al. (2006) have identified a weak absorption feature 
in the LMC sample at 5.8/im, attributed to carbonyl. This feature 
is not detected in our SMC sample. It could indicate an underabun- 
dance of CO with respect to stars in the LMC, as expected given 
that there is less O to start with. 
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4.3 Models of molecular opacity 

The top lines in Fig.|3]show the calculated model spectra of a num- 
ber of molecules. Similar plots in Zijlstra et al. (2006) are based on 
the line lists and models of J0rgensen et al. (2000). Newer line lists 
are now available for some simple molecules (Harris et al. 2002, 
2006), such as HCN and HNC. We also revised opacity curves for 
CO, C2 and CS from Zijlstra et al. (2006). 

The HCN/HNC linelist was computed using first principles 
quantum mechanics. It contains about 160000 rotation- vibration 
energy levels truncated at 18000 cm^^, and at J = 60. There are 
around 4 x 10* lines, but only 34 x lO'' of these are strong enough 
to contribute to opacity. Both the linear H-CN and H-NC geometric 
configurations are studied, as are some of the low lying delocalised 
states. The delocalised states are the states in which the H nucleus 
orbits the CN part of the molecule. The accuracy of the fundamen- 
tal vibrational transitions is about 3-4 cm~^; the error increases 
at higher energies. The transition intensities agree well with lab- 
oratory data. Incorporation of laboratory data into the linelist has 
significantly improved the accuracy of the frequencies of the tran- 
sitions between these low lying states. 

The opacity curves for CO, C2 and CS are calculated from 
absorption coefficients determined by Goorvitch (1994), Querci et 
al. (1994) and Chandra et al. (1995) respectively. 

Using this, we determined the opacity for the HCN, HNC, CO, 
CS and C2 as a function of wavelength. To determine the transmis- 
sion function, we then used the Lambert-Beer law (Banwell & Mc- 
Cash 1994), i.e. the solution to the ID equation of radiative transfer 
in the absence of an internal source. 

Fitting the spectral features to get column densities is be- 
yond the scope of this paper and will be done in a forthcoming 
paper. We can however determine which molecules are responsi- 
ble for the observed molecular bands. We thus calculated trans- 
mission curves using column densities for the different molecular 
species of lO^^cm"^ for the HCN, HNC, CS and C2 molecules and 
lO^'^cm"^ for CO. A temperature of 1750 K was used for all the 
models. The resulting transmittance curves are plotted at the top 
of Fig. (3) together with the transmittance curves from J0rgensen 
et al. (2000) for C2H2 and C3. This C2H2 model reproduces the 
broad photospheric absorption, but not the narrow component, aris- 
ing from a cool layer (J0rgensen et al. 2000). Cherchneff (2006) 
shows that acetylene reaches its peak abundance in the extended 
atmosphere, and its temperature is expected to be lower than the 
photospheric temperature (van Loon et al. 2006, Matsuura et al. 
2006). 



5 INDIVIDUAL STARS 
5.1 PMMR 52 

PMMR52 is a red supergiant of spectral type K (Prevot et al., 
1983). This is the only oxygen-rich star in the sample. The spec- 
trum is shown in Fig.|5] The observed spectrum of PMMR 52 shows 
the presence of spectral features of oxygen-rich dust (silicates at 
10/im and 18/^m). The 18/im silicate feature is strong. The lO-fim 
silicate band is remarkably weaker than observed in other O-rich 
stars. Normally, in O-rich stars, the 18-/im feature should decrease 
steeply after 19-20/im. But the one in PMMR 52 remains very flat. 

There are no signatures of crystalline silicates, which indicates 
that the silicates are largely in an amorphous phase. Some small ab- 
sorption features are observed at 6.2 and 7.5/^m. The 6.2/^m feature 
is due to water, while the one at T.Sfim might be attributed to ab- 
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Figure 8. The peculiar spectrum of GM 780. 



Table 3. Mid-infrared continuum bands for carbon stars for the so-called 
'Manchester system'. The last column gives the adopted flux con'esponding 
to zero magnitude. 



central A [/^m] 


A-range [fim] 


Fo [Jy] 


6.4 


6.25-6.55 


96.5 


9.3 


9.1-9.5 


45.7 


16.5 


16-17 


15.4 


21.5 


21-22 


9.1 



sorption by SiO, but this feature is very faint and might be due to 
noise. 



5.2 GM 780 

Fig. [S] shows the spectrum of GM 780. Compared to the other 
sources, C2H2 bands are noticeably absent. The weak band at 8 ^m 
is best fitted as CS, based on the models above. The 1 l-/im band is 
strong but its shape is significantly different from that of the other 
stars. There are no traces of silicates at 9.8 and 18/^m. GM780 is 
the brightest star at K in our sample; the bolometric luminosity is 
high for our sample but lower than those of the luminous IRAS- 
selected carbon stars (van Loon et al. 1999b). 

The J — K colour of this object, J ~K = 2.6, is quite red, and 
dust should be present around it. The dust temperature is relatively 
high, however, at 720 K. 

The object is discussed further below (Section [9]l where we 
argue the this star has a C/O ratio lower than that of other stars in 
the sample. 



6 COLOURS AND BAND STRENGTHS 

As discussed above, the spectra of the observed stars are dom- 
inated by molecular and dust bands. This makes the definition of 
the continua difficult. To avoid this and to determine colours, we 
use the Manchester system (Zijlstra et al. 2006, Sloan et al. 2006). 
This system, valid for carbon stars but not for oxygen-rich stars, 
defines four narrow bands, selected to avoid the molecular and dust 
bands, to determine the continuum. Table[3]shows the wavelengths 



Table 4. Photometry: colours measured using four narrow carbon stars con- 
tinuum bands. 



target 


r A~i rr\ o n 

[6.4]-[9.3] 


[16.5J — [21.5J 




[mag] 


[mag] 


NGC419LE 27 


0.067 ±0.018 




NGC419LE 18 


0.067 ±0.020 


0.126±0.098 


NGC419LE 35 


0.215 ±0.013 




NGC419LE 16 


0.325 ±0.021 


-0.052±0.062 


NGC 419 IRl 


0.496 ±0.011 


0.047±0.017 


NGC 419 MIRl 


1.343 ±0.012 


O.352±0.012 


IRAS 00554 


0.938 ±0.008 


0.232±0.021 


RAW 960 


0.199 ±0.022 


0.017±0.115 


ISO 00573 


0.255 ±0.017 


-0.012±0.029 


LEGC 105 


0.380 ±0.013 


0.034±0.069 


ISO 01019 


0.374 ±0.015 


-0.156±0.039 


ISO 00548 


0.399 ±0.013 


-0.084±0.025 


ISO 00549 


0.725 ±0.014 


0.199±0.012 


GM780 


0.388 ±0.016 


0.185±0.031 



used to defined the four narrow-bands and the adopted flux corre- 
sponding to zero-magnitude, determined to give a zero colour for 
a Rayleigh- Jeans tail. Using this system, we can determine contin- 
uum flux at 6.4, 9.3, 16.5 and 21.5/im, by integrating the observed 
spectra over the defined bands. The choice of the wavelengths used 
to define the continua is discussed in Zijlstra et al. (2006). Table|4] 
lists the measured continuum flux and colours for the observed car- 
bon stars. 

This method permits us to determine two colour temperatures, 
[6.4]-[9.3] and [16.5]-[21.5]. The [6.4]-[9.3] colour is an indi- 
cation of the optical depth, whereas the [16.5] — [21.5] colour can 
be use to determine the dust temperature. Table [6] lists the black- 
body temperatures derived from the [16.5] — [21.5] colours for the 
observed stars. For the bluest objects, this method gives lower 
limits lO^K) as the Rayleigh-Jeans limit is reached. This is 
the case for NGC 419 LE 16, NGC 419 IRl, RAW 960, LEGC 105, 
ISO 00573, ISO 01019 and ISO 00548. We thus do not give an es- 
timation of the dust temperature for these stars. Furthermore, the 
derived temperature can be affected by wavelength-dependent con- 
tinuum contamination from other sources in the beam. Thus, for 
the bluest stars, we can obtain temperatures that are not physically 
meaningful. 

Fig.|7]shows the [6.4]-[9.3] vs [16.5]-[21.5] colour-colour 
diagram obtained with the Manchester method. The stars from our 
sample form a well-defined sequence on this diagram. Some stars 
show a negative [16.5]-[21.5] colour (NGC 419 LE 16, ISO 00573, 
ISO 01019, ISO 00548). These stars are very blue and their spectra 
are typical of naked carbon stars. Two blue stars, NGC 419 LE 18 
and GM780, are offset in [16.5]-[21.5] colour by roughly 0.1 
mag. The spectrum of G1VI780 is very noisy in LL, so that its 
[16.5] — [21.5] colour has higher uncertainties than the other ob- 
jects. The good correlation observed between [6.4] — [9.3] and 
[16.5] — [21.5] indicates that the choice of our continuum wave- 
lengths is appropriate. 

For all the stars in our sample, we have measured the strength 
of the observed features. We define small wavelength ranges on 
the blue and red side of the features (Table [Sj, selected to avoid 
features, to define the continuum. We then fit line segments from 
both sides of the feature that defines the continuum. As the red 
edge of the MgS feature is outside of the IRS wavelength range. 
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Table 5. Wavelengths used to estimate the continua for the SiC and C2H2 
spectral features. 



Features 


A [fim] 


Blue continuum [fim] 


Red continuum [fj,m\ 


C2H2 


7.5 


6.08-6.77 


8.22-8.55 


SiC 


11.3 


9.50-10.10 


12.80-13.40 


C2H2 


13.7 


12.80-13.40 


14.10-14.70 



this method can not be applied to the MgS feature. The continuum 
under the MgS feature is thus assumed to be a blackbody with a 
temperature deduced from the [16.5] — [21.5] colour. After subtrac- 
tion of the continuum, we can measure the strength of each feature. 
For the molecular bands (C2H2 at 7.5 and 13.7/im), this strength 
is estimated as an equivalent width. For the dust emission features 
(SiC and MgS), we use a line-to-continuum ratio, defined as the in- 
tegrated flux of the band divided by the integrated underlying con- 
tinuum, over the wavelength range of the feature. The MgS feature 
extend beyond the edge of the spectral coverage: its strength relates 
only to the part blueward of 38/im, and the continuum is calculated 
as a black body (see Zijlstra et al. 2006). 

We also measured the central wavelength of the SiC band. 
This central wavelength is defined as the wavelength at which, af- 
ter continuum removal, the flux on the blue side equals the flux on 
the red side. Table |6] lists the measured strength and central wave- 
lengths of the observed features. The MgS feature central wave- 
length can not be measured from the available wavelength range. 
7 COLOUR-COLOUR DISCRIMINATIN BETWEEN C 
AND O-RICH AGB STARS 

Classification criteria for separating C-rich from 0-rich stars pro- 
posed by Egan et al. (2001) are based on the K— A versus J— K 
colours (where A is the MSX band centred at 8.3/im). These are 
now known not to work effectively (Zijlstra et al. 2006, Buchanan 
et al. 2006). Better separation is possible using different MSX 
bands (Whitelock et al. 2006) but in the Magellanic Clouds only 
the A-band was sensitive enough to detect the AGB stars. 

A variety of colour-colour plots based on IRAC and MIPS ob- 
servations, useful in distinguishing chemical type, are discussed in 
Blum et al. (2006). These, however, rely on the IRAC bands outside 
of the IRS wavelength coverage. We investigated whether a colour- 
colour diagram, based on Spitzer IRAC and MIPS colours, could 
discriminate 0-rich and C-rich stars using A > 5 fim. The equiv- 
alent IRAC and MIPS flux for the observed objects were derived 
from the IRS spectra and the transmission curves of the IRAC 5.8 
and 8.0^m filters and MIPS 24/im filter. As our sample contains 
just one 0-rich star, we added IRS spectra of evolved stars in the 
SMC and in the LMC from different programs. 

Fig|9] shows the distribution of the 0-rich stars (open sym- 
bols) and C-rich stars (black symbols) from the combined sample 
on a [8.0] — [24] vs [5.8] — [8.0] colour-colour diagram. The 0-rich 
and C-rich stars form two well-defined sequences on this diagram. 
The C-rich stars have redder [5.8] — [24] colour than 0-rich stars 
for the same [8.0] -[24] colour. FigfTOlexpl ains this effect. We have 
plotted the spectra of PMMR 52, an 0-rich star, and ISO 00549, a 
C-rich star, with the transmission of the IRAC 5. Sum (dashed line), 
8um (dotted line) and MIPS 24um (dotted-dashed line) bands su- 
perposed. The filter responses are affected by several molecular and 
dust bands. For 0-rich stars, the blue edge of the MIPS 24^m fil- 
ter overlaps with the 18/im silicate emission feature, and molecu- 
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Table 6. Strength of the molecular and dust features, in terms of either the equivalent width in microns, or the integrated line-to-continuum ratio (Sect.|6] The 
last column gives the continuum (black-body) temperature, derived from the [16.5] — [21.5] colour listed in Table|4] 



target 


EW (7.5 ^tm) 


EW(13.7Atm) 


L/C(SiC) 


Ac 


L/C (MgS) 


T(K) 


NGC 419 LE 27 


0.080 ±0.015 


0.057 ±0.009 


0.033±0.009 


11.96 ±0.17 






NGC419LE 18 


0.101 ±0.009 


0.151 ±0.038 


0.032±0.016 


11.46 ±0.20 




1038±430 


NGC 419 LE 35 


0.129 ±0.006 


0.079 ±0.006 


-0.009±0.004 


10.48 ±0.00 






NGC 419 LE 16 


0.158 ±0.011 


0.033 ±0.009 


0.056±0.009 


11.33 ±0.06 






NGC 419 IRl 


0.150 ±0.002 


0.044 ±0.005 


0.038±0.004 


11.05 ±0.07 






NGC 419 MIRl 


0.187 ±0.011 


0.066 ±0.005 


0.101 ±0.005 


11.38 ±0.05 


0.392 ±0.015 


409± 12 


IRAS 00554 


0.102 ±0.008 


0.079 ±0.004 


0.119±0.008 


11.28 ±0.08 


0.276 ±0.026 


589±44 


RAW 960 


0.230 ±0.009 


0.051 ±0.010 


0.057±0.010 


11.20 ±0.11 






ISO 00573 


0.148 ±0.005 


0.078 ±0.013 


0.056±0.011 


11.19 ±0.18 






LEGC 105 


0.116 ±0.003 


0.062 ±0.006 


0.019±0.004 


11.01 ±0.11 






ISO 01019 


0.151 ±0.005 


0.032 ±0.013 


0.033±0.009 


11.26 ±0.13 






ISO 00548 


0.138 ±0.003 


0.068 ±0.002 


0.068±0.004 


11.27 ±0.04 






ISO 00549 


0.047 ±0.013 


-0.041±0.010 


0.061±0.012 


11.21 ±0.30 




674±34 


GM 780 


0.149 ±0.010 


0.006 ±0.007 


0.289±0.011 


11.13 ±0.05 




724± 94 




Figure 9. [8]— [24] vs [5.8] — [8.0] colour-colour diagram of AGB stars in 
the LMC and SMC. Black symbols represent AGB C-rich stars. Open sym- 
bols represent AGB O-rich stars. Squares coiTespond to stars in the SMC 
and triangles to stars in the LMC. 



Figure 10. Spectra of two stars with the same [8] — [24] colours. No flux 
scaUng has been applied. ISO 00549 (top) is a C-rich star. PMMR52 (bot- 
tom) is an O-rich star. The transmission of the IRAC 5.8/jm (dashed line), 
8fim (dotted line) and MIPS 24/^m (dot-dashed line) filters are overlayed. 



lar absorption by SiO affects the 8/im band. For carbon stars, the 
5.8/im band is suppressed by molecular absorptions. (We note that 
this IRAC band extends beyond the blue edge of our spectra.) A C- 
rich star with the same [8] — [24] colour as an O-rich star will have 
a continuum emission redder than the O-rich star. Such a diagram 
has been made by Buchanan et al. (2006) using their LIVIC sample 
containing redder stars. 

Fig[T0]thus indicates that a [8]-[24] vs [5.8]-[8.0] diagram 
can be used to discriminate between C and O-rich stars with colours 
within the range of the objects present in our sample ([8] — [24] <4 
and [5.8] — [8.0] < 1.4). Furthermore, the sample used in Fig|9]con- 
tains stars in the LIVIC (triangles) and the SIVIC (squares). The SIVIC 
and LMC follow the same sequences, indicating that the difference 
of location of C-rich and O-rich stars in the colour-colour diagram 
is independent of metallicity. 



8 AVERAGE SPECTRA FOR SMC AND LIVIC STARS 

Fig. [TT| shows the averaged spectra of the red SIVIC sources (top, 
solid line) and blue SIVIC sources (bottom, solid line), where 'red' 
and 'blue' are defined as in Zijlstra et al. (2006). To compare these 
spectra with LIVIC spectra, we overplot the LIVIC spectra on this 
figure (dotted lines). Those LIVIC spectra are taken from Zijlstra et 
al., and are respectively the spectra with weak SiC and no IVIgS in 
Fig. 19 in that paper. The right panel shows the ratio of the SIVIC 
averaged spectra and the LIVIC averaged spectra. 

The LIVIC and the SIVIC spectra show the same features, but 
the differences are significant especially for the 'red' sources. The 
differences are associated with the C2H2 and the SiC bands. The 
fact that the acetylene features show as absorption in the ratio spec- 
tra indicates that these bands are relatively stronger in the SIVIC. 
The 13.7/im is clearest seen, indicative of the colder circumstellar 
molecules. The 11.3/im band also shows in absorption, but as this 
is an emission feature, the negative footprint here shows it to be 
weaker in the SIVIC. 

The 'blue' sources do not show easily interpretable differ- 
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Figure 11. On the left, averaged spectra from our SMC sample (solid lines). The upper solid line shows the averaged spectrum of the red sources. The lower 
solid line shows the averaged spectrum of the blue sources. The dotted lines show the equivalent averaged spectra from the Zijlstra et al. LMC sample. The 
right panel shows the SMC averaged spectrum divided by the LMC averaged spectrum for the red sources (top) and the blue sources (bottom). 



ences. The broad lO^m absorption may contain a blue contribution 
from a molecular band (e.g. Zijlstra et al. 2006), in addition to SiC. 
The redward continuum may be affected by the photospheric broad 
14/im C2H2 band. 



9 SILICON-CARBIDE DUST 
9.1 Strength 

Figll II clearly indicates that the 11.3 /xm SiC feature is stronger 
in the LMC than in the SMC, even if this feature is rather weak 
in blue sources. The divided spectra do not show any difference in 
the shape of the SiC feature in the SMC and the LMC, but only in 
the band strength. Thus, whereas the emission strength varies, the 
emission shape is constant. The composition of the SiC grains is 
likely the same for the two environments. 

Fig[T2] shows the strength of the SiC feature as a function of 
the [6.4] — [9.3] colour. To study the effect of metallicity on the 
strength of the SiC features, we overlaid the measured strength of 
this feature on the Galactic and SMC sample of Sloan et al. (2006) 
as well as our LMC sample (Zijlstra et al. 2006). Sloan et al. (2006) 
find weaker SiC and MgS features in the SMC than in the Galaxy, 
related to the lower abundance of Si, S and Mg in the SMC. 



Our SMC sample confirms that the strength of the SiC feature 
is lower in the SMC than in the LMC and the Galaxy. However, 
there are two noticeable exceptions from the Sloan et al. sample 
(filled triangles). 

The distribution of strength with colour (optical depth) is fun- 
damentally different in the Galaxy compared to the SMC. In the 
SMC, there is a clear trend of increasing SiC strength with optical 
depth (left panel of Fig|12t. In the Galaxy, SiC increases rapidly 
up to [6.4] — [9.3]^ 0.5, followed by a decline towards the SMC 
relation. The difference in the Galactic and Magellanic Cloud dis- 
tributions is very clear in the right panel of Fig. [T2] Note that one 
SMC star follows the Galactic sample: GM780, located outside 
of the scale of the left panel (and therefore not shown). This object 
(see Section [5!2l l has a strong SiC feature but with an unusual shape 
and rather weak dust continuum. We return to it below. 

The rapid increase of SiC L/C for Galactic stars at very low 
[6.4] — [9.3] colour (Fig. I12t suggests that in the Galaxy, the SiC 
feature is present while the dust excess is still low. The subsequent 
decline in L/C ratio may include effects of optical depth within 
the feature (Speck et al. 2005), or it may be due to the increasing 
dust continuum. The LMC stars increase much more slowly with 
optical depth, and reach the peak around [6.4] — [9.3] « 1.0; red- 
ward they are close to or a little below the Galactic stars. For the 
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Figure 12. The strength of the SiC feature as a function of the [6.4] — [9.3] colour in our SMC sample (left panel). The right panel combines this with stars 
from Sloan et al. SMC and Galactic samples: filled squares represent our SMC sample, tiiangles Sloan et al. (SMC), open circles Sloan et al. (Galactic) sample 
and filled circles Zijlstra et al. (LMC). 



SMC stars, there is no well-defined peak but most SMC stars with 
[6.4]-[9.3]< 1.0 have little or no SiC. 



9.2 Condensation sequence 

The difference between the SMC and the Galaxy can be understood 
if for Galactic stars, SiC condenses first (i.e. at higher tempera- 
tures) and carbon dust later. For the LMC stars, the two form close 
together, while for the SMC stars, SiC forms last. Thus, there is a 
reversal in condensation sequence over this range in metallicity. 

The condensation temperature of graphite is around 1600K 
whilst for SiC Tc ~ 1400 K. However, these are composition and 
density dependent. For high partial pressure, condensation occurs 
at higher temperatures. The effect is shown by Bernatowicz et al. 
(2005). For increasing C/0 ratio, graphite condenses at higher T. 
With decreasing Si abundance, Tc (SiC) will decrease. The effect 
is that SiC condenses first at high density, low C/O ratio (< 1.1) 
and high Si abundance. Fig. [T2] suggests these conditions are met 
for Galactic stars. For LMC stars, SiC condensation is delayed, 
while for SMC stars in many cases SiC does not condense at all. 
The increasing C/O ratio also aids this process. Chigai & Ya- 
mamoto (1993) also suggest that under different conditions, SiC 
and graphite may form together, only SiC grains may form or only 
graphite forms. 

Following this, we interpret GM 780 as a case where the den- 
sity in the dust-formation region is high and the C/O ratio low, lead- 
ing to SiC only. The relatively blue Manchester colour of this object 
is due to the resulting lack of dust. The lack of clear C2H2 features 
is consistent with this interpretation. 

If the condensation sequences discussed above are correct, the 
SiC abundance may be limited by Si in the Galaxy (where SiC con- 
denses first) but by the C/O ratio in the Magellanic Clouds (where 
much of the carbon goes into the dust before SiC can form). 



9.3 Central wavelength 

Fig.[T3]shows the apparent central wavelength of the SiC feature as 
a function of the [6.4] — [9.3] colour. The right panel displays our 
measurements together with the other samples mentioned above. 
Our sample shows that for the stars with [6.4] — [9.3]<0.5, there 
is a clear trend of decreasing central wavelength with increasing 
colour, while the opposite is observed for redder stars. The dashed 
line shows the fit proposed for Galactic stars (Sloan et al. 2006). 

This shift could be due to two factors (Zijlstra et al. 2006). The 
broad C2H2 absorption band centred at 13.7/im is close to the SiC 
feature and the red continuum wavelengths we use for SiC neces- 
sarily fall within this band. Thus stronger C2H2 absorption shifts 
the SiC central wavelength to the red. Speck et al. (2005) also sug- 
gested that self-absorption of SiC can also induce a shift toward the 
red of the SiC feature. 

The shape of the SiC feature appears to be the same for all 
stars (e.g. Zijlstra et al. 2006). This argues against wavelength- 
dependent optical depth effects (Speck et al. 2005) for our samples. 
We note that the stars studied here have relatively low optical depth 
compared to stars with known SiC absorption. 



10 MAGNESIUM SULFIDE 

Longward of 15/im, only the MgS feature is seen. The shape of this 
feature is similar in the SMC and the LMC. Furthermore, we have 
shown (Sect l4.lt that this feature is present in the envelopes of the 
reddest stars, i.e. the stars with the lowest dust temperature. This 
is explained by the formation process of MgS which starts around 
600K and is complete around 300K. 

The feature is only seen for two stars in the SMC sample, but 
this appears to be due to the dust temperatures. Fig. [14] shows the 
comparison between the SMC and the LMC. There is little dif- 
ference between them, and it appears the ratio between MgS and 
carbon dust is less sensitive to the metallicity than is the SiC versus 
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carbon dust. This may indicate that the MgS condensation is limited 
by the available dust surface rather than the element abundances. 

The shape of the continuum subtracted MgS feature is shown 
in Fig[T5] The MgS band of NGC 419 MIRl has a relatively flat 
shape, while the band observed in IRAS 00554—7351 diminishes 
toward the red. These shapes are similar to the shapes of the MgS 
features observed in the LMC (Zijlstra et al. 2006). The study of 
this feature in Galactic sources has shown that it could be resolved 
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Figure 15. Shape of the MgS feature, continuum subtracted, for SMC 
carbon-rich AGB stars. 



in two subpeaks centred at 26 and 33/im (Volk 2002). Those sub- 
peaks are not observed in our sample, but this is not conclusive due 
to the low signal. 



11 GAS PROPERTIES 
11.1 C2H2 

Fig. 1161 shows the equivalent width of the C2H2 band at 7.5/im, 
for SMC, LMC and Galactic stars. Due to a lower dust-to-gas ra- 
tio at lower metallicity, stars with the same gas mass-loss rate will 
have a lower optical depth in the SMC than in the LMC and the 
Galaxy. If we use the [6.4] — [9.3] colour as a tracer of the dust 
mass-loss rate, the strength of the 7.5/im band is comparable in the 
Galactic sample and our SMC sample. However, the [6.4] — [9.3] 
colour will be different for a given gas mass-loss rate. Zijlstra et 
al. (2006) show that the effect of metallicity on the optical depth 
(lower metallicity stars have lower optical depth) is stronger for 
redder stars (see Fig. 10 in their paper). Most of the SMC stars in 
our sample have [6.4] — [9.3]<0.4. To be compared with stars with 
similar gas mass-loss rates, we should, according to the previous 
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Figure 17. The strength of the C2H2 7.5/xm feature as a function of the 
Sic feature. Squares represent SMC stars from our sample and Sloan et al. 
sample, black circles Zijlstra et al. LMC sample and open circles Sloan et 
al. Galactic sample. 

work, compare those stars with stars with [6.4] — [9.3]« 0.5 and 
0.6 in the LMC and the Galaxy, respectively. 

By applying this shift, we observe that the T.Sfim feature is 
stronger in the SMC than in the LMC and the Galaxy. This feature 
appears to be stronger in lower metallicity environments. This has 
previously been found for the LMC (Zijlstra et al. 2006, Matsuura 
et al. 2006) and the SMC (Sloan et al. 2006). The stars in our sam- 
ple are bluer than the previous observed stars, and we confirm that 
the C2H2 feature is stronger in low metallicity environments even 
when the stars are almost naked stars. 

As noticed before in the LMC and SMC, this effect is more 
clearly seen for the 7.5/im feature than for the 13.7/im one. 

The divided spectra in Fig.[TT]present absorption bands at 7.4 
and 13.7/im. This confirms that the C2H2 features are stronger in 
the SMC than in the LMC. 

Fig[T7] shows a trend of increasing C2H2 strength with in- 
creasing SiC strength in the SMC. Such a trend is also observed 
in the LMC but the slope is different. Whether this trend is present 
for the Galaxy is less obvious. The SiC condensation sequence dis- 
cussed above is consistent with this: in the Galaxy, SiC is limited 
by Si, and in the Magellanic Clouds by the C/0 ratio. This predicts 
a correlation with C2H2 only for the Magellanic Clouds. 

The distribution in Fig[T7] clearly shows a separation into 
SMC-LMC~Galaxy, from left to right. The ratio of C2H2 (gas) 
over SiC (dust) is higher in the SMC, then in the LMC and the 
Galaxy. This shows an increasing gas-to-SiC dust mass ratio with 
decreasing metallicity. 

11.2 CO 

In all the spectra, the flux drops sharply shortward of 6/im. This is 
due to absorption by the CO molecule (j0rgensen 2000), C2 and/or 
C3 (Zijlstra et al. 2006), as shown by the spectra at the top of Fig. 
[S] The slope of this feature is steeper in the LMC than in the SMC. 
We attribute this to an underabundance of CO in the SMC with 



respect to the LMC. For carbon stars, the CO abundance is limited 
by oxygen. Oxygen is a measure for the original metallicity: the 
CO band is therefore expected to vary with metallicity. Note that 
the responsivity at the blue edge of the spectra is low. 



12 CARBON VERSUS OXYGEN-RICH STARS 

A noticeable finding for both the LMC and the SMC is that the 
Spitzer targets are found to be strongly biassed towards carbon 
stars. The selection criteria were relatively insensitive to the chem- 
ical type, and the bias is therefore intrinsic rather than a selection 
effect. Apart from stars observed as a result of acquisition errors, 
we found 14 carbon stars and no oxygen-rich stars. 

Sloan et al. (2006) use 2MASS colours as part of their se- 
lection criteria. This gave a sample of 22 carbon stars versus 
10 oxygen-rich stars (excluding some other sources). The use of 
2MASS requires low circumstellar reddening and favours stars with 
silicate dust where the reddening per unit dust mass is lower. 

The LMC sample of Zijlstra et al. (2006) shows 28 carbon 
stars versus one oxygen-rich star (excluding one high-mass proto- 
star). Buchanan et al. (2006) observed high-luminosity stars and 
found 16 C-rich versus 4 0-rich AGB stars. Thus, the dominance of 
carbon stars is confirmed in all surveys of mass-losing stars, but the 
ratio depends on the characteristics of the survey. The proportion of 
C stars is lower in high luminosity, lower-optical depth surveys. 

Cioni & Habing (2003) show that for optically visible evolved 
stars in the Magellanic Clouds, the ratio of C over M stars is ap- 
proximately 3 for both clouds. Blum et al. (2006) surveying the 
full population of AGB stars in the LMC, identify 17500 O-rich 
and 7000 C-rich stars. The higher ratio of O to C stars found by 
Blum et al. may suggest that they include K-type stars, earlier on 
the AGB. Cioni & Habing (2003) estimate that there are about 2250 
C stars in the SMC and 10000 in the LMC, the latter being roughly 
consistent with the Blum et al. nuinbers. 

The high ratio (~3) of C to O-rich stars is consistent with 
the expectation that in the Magellanic Clouds, most stars become 
carbon stars early on the AGB, before the onset of the high mass- 
loss rate phase. There are almost no O-rich stars in the present and 
the Zijlstra et al. (2006) samples of lower luminosity, high mass 
loss rate, AGB stars in the Magellanic Clouds. Only among the 
highest mass, higher luminosity stars on the AGB do O-rich stars 
contribute significantly to the population of mass-losing stars." 



12.1 Carbon star evolution 

Fig. [6] shows that in the LMC, the carbon star luminosity function 
exhibits a flat peak with a width of 0.3 mag. This traces the part of 
the AGB evolution where essentially all stars have become carbon 
stars: the flat distribution results from the linear change of mag- 
nitude with time. (Stars more massive than ~ 4 M0 may avoid a 
carbon star phase due to hot bottom burning, but there are few such 
stars.) At the high-luminosity end, the number of known stars drops 
because of increasing obscuration and stars leaving the AGB. 

In the SMC, the width of the flat peak is about 0.6 mag. Over 
10'' yr, a star brightens by approximately 1 mag on the AGB. These 
numbers therefore predict that the optical carbon star phase lasts 
3 X lO'^ to 6 X 10^ yr, for the LMC and the SMC respectively. 
These are lower limits to the total life time of the carbon star phase 
on the AGB. 

The SMC carbon star distribution peaks at lower luminosity 
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Figure 16. The equivalent width of the C2H2 band at 7.5/ini (left) and 13.7/im (right) as a function of the [6.4] — [9.3] colour. Filled squares represents our 
SMC sample, triangles the Sloan et al. SMC sample, open circles the Sloan et al. Galactic sample and filled circles the Zijlstra et al. LMC sample. 



than in the LMC. The longer life time in the SMC is therefore due to 
the stars becoming carbon-rich earlier, after fewer thermal pulses. 

The Spitzer targets all are at the high-luminosity end, and trace 
the final stage of their evolution. The Manchester colours of the 
SMC stars indicates mass-loss rates of 10~^ M© yr~^ or higher 
(Zijlstra et al. 2006, their Fig. 10). Assuming that a typical carbon 
star expels 0.2 Mq during the superwind, indicates a duration of 
this phase of order ~ 10* yr. The combined SMC samples have 
approximately 30 superwind stars, or roughly 1 per cent of the to- 
tal number of carbon stars predicted by Cioni & Habing (2003). 
This suggests a life time of the carbon star phase of ^ 10*' yr. 
For comparison, the thermal-pulsing AGB lasts 1-2 x lO*" yr for 1- 
3 M0 stars. Thus, SMC stars become carbon stars during the first 
half of the thermal-pulsing phase. 



This leads to the suggestion that at low metallicity, stellar mass 
loss is mostly carbon-rich, as at the luminosity where stellar pulsa- 
tions develop an extended atmosphere leading to molecule and dust 
formation, the stars have already become carbon rich. Only in old 
stellar systems (globular clusters. Galactic halo), where the stellar 
masses are insufficient for third dredge-up, is an oxygen-rich wind 
expected. 

The lack of silicon implies that the dust input into the ISM 
from LIMS in low metalUcity systems would be lacking in sili- 
cates and silicate-carbides, but be dominated by carbonaceous dust, 
a 'soft' (mineral-poor) type of dust (Zijlstra 2006, Matsuura et al. 
2006, Groenewegen et al. 2007). 



12.2 Mass loss and metallicity 

A clear result from the Spitzer surveys is the enhancement of C2H2 
in low metalUcity environments. This has previously been foimd 
from ground-based spectra. Van Loon et al. (1999b) and Matsuura 
et al. (2002, 2005) argue that this is related to a higher C/0 ra- 
tio in these stars. The C/O ratio is enhanced by two effects: (i) a 
lower oxygen abundance means that less carbon is locked up in CO 
(e.g. Lattanzio & Wood 2004), and (ii) lower metallicity promotes 
more efficient third dredge-up of carbon (Wood 1981, Vassiliadis & 
Wood 1993) with stars experiencing more thermal pulses (Lawlor 
& MacDonald 2006). 

This seems to be confirmed by the finding that LMC car- 
bon stars have higher abundances of C2H2 than Galactic stars. In 
contrast, oxygen-rich stars show low molecular abundances at low 
metallicity (Zijlstra 2006). 

Spitzer shows that in the LMC and SMC, high mass-loss stars 
are heavily dominated by carbon stars, unlike in the Galaxy (this 
paper, Zijlstra et al. 2006, Sloan et al. 2006, van Loon et al. 2006, 
Buchanan et al. 2006, Blum et al. 2006). In contrast, the sole 
oxygen-rich AGB star in the sample of Zijlstra et al. has a rather 
moderate mass loss. Only the most Iimiinous LMC stars show an 
0-rich population with thick dust shells. 



13 CONCLUSION 

We have presented a Spitzer spectroscopic survey of 14 carbon-rich 
AGB stars in the SMC. The bolometric magnitudes of the observed 
stars indicates that their initial masses are in the range \-AMq. We 
also presented spectra of an accidentally-observed K supergiant and 
an HII region. 

Our spectra covers the range 5-38 /im. Molecular bands due to 
C2H2 at 7.5 and 13.7um are observed in most of the C-rich stars. A 
weak absorption band at 14.3 /um, attributed to HCN, is observed 
in the two reddest stars of our sample. The 5.8-micron band of car- 
bonyl observed in the LMC is not observed in our SMC sample. 
This could be due to an underabundance of CO in the SMC with 
respect to the LMC. The acetylene bands provide clear evidence 
for a higher C/O ratio in low metallicity carbon stars. 

Using the "Manchester System", we determined the contin- 
uum flux for the observed stars using four narrowbands centred at 
6.4, 9.3, 16.5 and 21.5/im. This enables us to determine two colour 
temperatures, [6.4]— [9.3] and [16.5]— [21.5]. These colours are in- 
dicators of the optical depth and dust temperatures respectively. 

We find evidence for a different behaviour of the SiC feature 
with metallicity: in the Galaxy, it appears for stars with little dust 
excess, while for the SMC it only appears for stars with cool dust. 
The LMC is intermediate. We suggest that this is due to a differ- 
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ent formation sequence, where at solar metallicity SiC precedes 
graphite, while at SMC metallicity and its corresponding higher 
C/0 ratio, graphite forms first. 

The MgS feature at ~30 /im has strength comparable to LMC 
stars, suggesting its abundance relative to that of amorphous carbon 
is less metallicity dependent. We suggest its formation is limited 
by the available surface area of pre-existing dust grains, rather than 
elemental abundances. The shape of the MgS and SiC features are 
similar in both galaxies, but the SiC feature is much stronger in the 
LMC than in the SMC. 

We show that a colour-colour diagram using Spitzer IRAC 
and MIPS filters ([5.8]~[8.0] vs [8.0]-[24]) is able to separate 
O-rich and carbon stars. This had been shown by Buchanan et al. 
2006 for redder stars; we confirm that it is also valid for blue stars 
([8] — [24])<3). Furthermore, we show that this separation is inde- 
pendant of metallicity. This suggests that it is possible to discrim- 
inate C and O-rich AGB stars in other galaxies using only Spitzer 
photometry. 

Carbon stars dominate the population of high mass loss rate 
AGB stars in the SMC. We derive a life time of the carbon star 
phase of 5J 6 X 10'' yr, with the superwind phase lasting 10* yr. 
Thus, stars become carbon-rich early on the thermal-pulsing AGB, 
long before the onset of the superwind. The LMC stars become 
carbon-rich later in their evolution, and last as carbon stars for ^ 
3 X 10^ yr. 
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